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The conceptual and preliminary design processes for aircraft with large shape changes 
are generally difficult and time-consuming, and the processes are often customized for a 
specific shape change concept to streamline the vehicle design effort. Accordingly, several 
existing reports show excellent results of assessing a particular shape change concept or 
perturbations of a concept. The goal of the current effort was to develop a multidisciplinary 
analysis tool and process that would enable an aircraft designer to assess several very 
different morphing concepts early in the design phase and yet obtain second-order 
performance results so that design decisions can be made with better confidence. The 
approach uses an efficient parametric model formulation that allows automatic model 
generation for systems undergoing radical shape changes as a function of aerodynamic 
parameters, geometry parameters, and shape change parameters. In contrast to other more 
self-contained approaches, the approach utilizes off-the-shelf analysis modules to reduce 
development time and to make it accessible to many users. Because the analysis is loosely 
coupled, discipline modules like a multibody code can be easily swapped for other modules 
with similar capabilities. One of the advantages of this loosely coupled system is the ability to 
use the medium- to high-fidelity tools early in the design stages when the information can 
significantly influence and improve overall vehicle design. Data transfer among the analysis 
modules are based on an accurate and automated general purpose data transfer tool. In 
general, setup time for the integrated system presented in this paper is 2-4 days for simple 
shape change concepts and 1-2 weeks for more mechanically complicated concepts. Some of 
the key elements briefly described in the paper include parametric model development, 
aerodynamic database generation, multibody analysis, and the required software modules as 
well as examples for a telescoping wing, a folding wing, and a bat-like wing. 


I. Introduction 

A ll aircraft undergo some degree of shape change to accomplish their mission. This is commonly achieved by 
moving control surfaces. However, several aircraft have exhibited larger than normal shape changes that 
include variable sweep, variable span, and/or folding wing concepts. Some of the early research aircraft that 
exhibited large-shape change included the X-5, XB-70 Valkyrie, and the Oblique wing. Operational aircraft that 
have had the ability to change their shape significantly during their mission include the B-l Lancer, F-8, F-lll, F- 
14, Mirage G, MiG-23, MiG-27, Sukhoi Su-7, and Tu-22. In each of these examples, the specified mission was more 
diverse or broader than conventional design, driving the wing design to a highly reconfigurable or large-shape 
change concept. Large-shape change concepts nearly always have associated design “penalties” such as added 
weight or complexity. However, when overall system performance and mission requirements are assessed, large- 
shape change concepts can be a viable approach for some missions - particularly missions that may be described as 
combinations of several missions with diverse speed, altitude, or take-off and landing requirements. 

Over the past decade, several research efforts have focused on developing or assessing a new breed of multi-mission 
aircraft and related innovative technologies that include large shape change 1 . Technologies that can address multi- 
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mission requirements range 2 from little or no shape-change, such as micro flow control and variable bumps 3 to 
moderate shape change such as variable airfoil shape (including camber change 4 ) and wing twist 5 ; to large scale 
shape 6 change such as variable span 7 , sweep 8 , and fold angle 9 mentioned in the examples above. In the literature, 
this suite of technologies has many names including: reconfigurable, multifunctional, adaptive, or morphing 
technologies. This report will refer to advanced technologies that may address multi-mission requirements as 
“morphing.” Morphing may also be defined as efficient, multi-point adaptability. Technological advances and 
innovation in several fields, such as materials, structures, aeroelasticity, flow control, aerodynamics, and flight 
control, have enabled many of these novel morphing concepts to come closer to reality for aircraft in the 21 st 
century. 

With potential new capabilities in hand, today’s aircraft designer has new opportunities as well as new 
challenges, for the promises of morphing benefits are not without a penalty. Shape change typically requires 
additional hardware, weight, and/or complexity. One major weight penalty is the additional actuator weight, which 
is a function of the mechanisms used and the energy required to morph. Another design challenge is determining 
which technology or technologies to use for shape change, whether it is for little or no shape change or large shape 
change. Designers of morphing vehicles must also optimize the use of the technologies or suite of technologies and 
conduct trade-off studies to weigh benefits and penalties. The system benefits must significantly outweigh the 
penalties for the new morphing concept to survive the conceptual design phase. Further, assessment of shape-change 
concepts requires analysis of the air vehicle at different flight conditions and shapes, and each condition may require 
a different analysis model (e.g., aerodynamic mesh). Analysis models for small and moderate shape changes can be 
generated with slight modifications of the un-morphed analysis models. However, analysis models for large-scale 
shape changes may need to be regenerated for every morphed stage. This time-consuming model regeneration 
process often precludes conducting multiple analyses early in the design stage (during conceptual design) unless the 
process is automated for different morphing stages. Yet, conceptual design phase is where the aircraft designer has 
greatest flexibility to trade a variety of shape-change concepts. Unfortunately, due to the complexities mentioned 
previously, many conceptual design efforts using morphing concepts tend to focus on perturbations about one 
particular type of shape change concept (e.g. telescoping) and not on comparing several different concepts. 

The research effort reported herein focuses on developing a capability to assess several diverse morphing 
concepts (focusing on large shape change) with greater fidelity and efficiency. This approach essentially blends 
some of the tools normally associated with preliminary design (such as second-order performance analysis) with the 
efficiency and flexibility normally associated with conceptual design. This was achieved by integrating several 
existing analysis capabilities, which addressed some of the common challenges with multi-disciplinary analysis 
(MDA) for morphing aircraft. In bringing these analysis capabilities together, the authors found that the initial set-up 
was the most time consuming part of the process because we had to understand the assumptions and limitations of 
the different analytical tools. Once this was accomplished, the exercise of assessing a new morphing concept and 
getting second order performance results (e.g., loads, roll, and roll rate) was extremely quick - typically only a few 
days. This process would normally take months using the conventional preliminary design approach. 

In report, the overall MDA design approach is presented as well as a detailed description of each component of 
the design: model generation, data transfer, aerodynamic database generation, and multibody analysis. To evaluate 
the integrated design approach described herein, three generic and diverse shape change configurations were 
analyzed, and the results are provided for reference. 

II. Multidisciplinary Analysis (MDA) for Aircraft with Large-Shape Changes 

Design and assessment of reconfigurable aircraft introduce several challenges in the conceptual and preliminary 
design phases. Even for more traditional design approaches, stepping through the process of developing the simple 
sketches that initiate a conceptual design and completing the detailed analysis of a preliminary design can be very 
costly and time-consuming. Though preliminary design typically results in a better understanding of the benefits and 
penalties of a particular concept design, at this point in the design process, it is often too late to make significant 
changes to the design. During the NASA Morphing Program and the DARPA Morphing Aircraft Structures (MAS) 
program, a common question was frequently posed: How can we quickly assess several different morphing 

concepts without having to spend considerable resources and time doing a detailed design for each concept? Several 
challenges were identified and are described below. The current research effort was focused on addressing some of 
the challenges to enable a user-friendly, medium-fidelity, and yet time-efficient approach to MDA for aircraft with 
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large shape changes. This approach is centered on automating time-consuming steps, using existing off-the-shelf 
software modules, and integrating these modules in a loosely coupled manner. This way aircraft designers can use 
their existing tools along with this integrated tool for a rapid assessment of the types of and the amount of shape 
change that would address multi-mission requirements. 

One of the classic challenges of MDA is increasing fidelity while maintaining or increasing the number of 
disciplines and addressing the efficient and accurate data transfer among disciplines. Typically, it is very time- 
consuming and can be tedious to address all of these requirements because multiple models from different 
disciplines must be generated that include an appropriate level of physics and geometry characteristics. In addition, 
multi-point optimization may require multiple models for the aerodynamic analysis and structures. Morphing 
concepts may also require additional non-traditional disciplines such as multibody analysis (included in this study) 
or kinematics, etc. Thus, model generation for MDA of morphing concepts can be a significant challenge. 
Correspondingly, accurately moving data between the different disciplines is complex and often done in an ad-hoc 
manner, with no consideration for constraints such as conservation of forces, moment, and virtual work. An 
additional complication for vehicles with large shape change is that the uniqueness of the configuration often drives 
MDA analysts to develop a special-purpose tool (e.g. finite element analysis) to analyze the different shape changes. 
This special-purpose tool then uses a customized geometry description of a specific morphing concept. While this 
approach usually yields excellent results, analyzing a different concept requires a substantial change to the system. 
In addition, this leads to comparing different morphing concepts with different analysis systems, introducing 
inconsistency in the results. To address these challenges several steps were taken in the current study: 


1. We have used a parameterized model with a few critical parameters describing the geometry, flight 
conditions, and the morphing concepts. This parameterization allows one to analyze a diverse set of shape 
change concepts and enables automation of model generation and analysis. In the examples reported herein, 
an automated aerodynamic analysis was developed and integrated. In addition, automated finite element 
model generation and analysis were developed but not integrated. 

2. We used a general purpose and validated data transfer tool that is accurate and fully automated. This 
capability simplified the integration process. 

3. We used a loosely coupled integration method, which enables interchangeability of software modules with 
similar capabilities and allows the use of up-to-date analysis modules. In addition, the software modules are 
encapsulated and there is no need to modify them. This also enables the use of off-the-shelf software and 
the inclusion of more disciplines. In our examples, we were able to compare the results of two different 
multibody codes. 

4. We used off-the-shelf software modules to accelerate our development process. This put an additional 
burden on the initial integration effort. However, we believe this burden is a worthwhile effort as we were 
able to avoid writing large pieces of software and gained flexibility in being able to include different codes. 


The overall integrated 
approach is capable of 
analyzing a diverse set 
of morphing concepts, 
hence eliminating any 
inconsistency due 
dissimilar and 

disparate systems. 
Figure 1 shows the 
current integrated 
approach for design 
and analysis of 
morphing concepts. 
The approach includes 
geometry generation, 
aerodynamic analysis, 
structural analysis, 
multibody analysis, 
and data transfer 



Fig. 1 Overall design process. 
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among various disciplines. Modules for structural model generation and analysis have been developed, but they have 
not been integrated yet. The commercial software MATLAB® (MathWorks^) is used to generate the aerodynamic 
meshes; Discrete Data Transfer Between Dissimilar Meshes (DDTBDM 10 ) is used to transfer data among 
aerodynamics, structures, and the multibody dynamics models; Virtual Lab® (LMS^) and DYMORE 11 are used for 
dynamics analysis; and CMARC® (AeroLogic §§ ) is used for aerodynamics analysis. This is a loosely coupled 
integrated analysis, and one can easily swap any of these modules with other software modules with similar 
capabilities. 


III. Parametric Model Generation 

The use of parametric models is an essential element of our proposed approach. It allows for simplification and 
automation of the model generation for a diverse set of morphing concepts. The parametric model is based on a 
wing, which is decomposed into several panels. Each panel is parameterized as shown in Fig. 2 with five variables: 
1) root chord C r , tip chord C t , aspect ratio AR, 
quarter-chord sweep angle X, and wing fold 
angle 0. With this set of parameters, it is 
possible to generate a large number of 
morphing wings including the three examples 
shown in Fig. 3. Tables 1-3 show the 
parameters used in our examples and Figs. 4-5 
show the wing planforms. Table 4 shows the 
mass and inertia properties for the telescoping, 
folding, and bat-like wings, and these wings 
sized to be consistent with empirical vehicle 
sizes. 


c* 



Fig. 2 Geometry parameterization for one panel. 



a) Telescoping wing 



b) Folding wing 



c) Bat-like wing 

Fig. 3 Sample morphing wings (left: un-morphed, middle: 50% morphed, right: fully-morphed). 


^ MATLAB is a trademark of Math Works (http://www.mathworks.com) 

** Virtual Lab is a trademark of LMS International (http://www.lmsintl.com) 
§§ CMARC is a trademark of AeroLogic (http://www.aerologic.com) 
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Fig. 4 Planform view of telescoping and folding wings. 
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Fig. 5 Planform view of bat-like wing. 


Telescoping Wing 


Initial 

Final | 

Panel # 

1 

2 

1 

2 

Cr(meter) 

4 

1 

4 

1 

Ct(meter) 

1 

1 

1 

1 

AR 

1.2 

0 

1.2 

3 

X (deg) 

37 

0 

37 

0 

Fold Angle (deg) 

0 

0 

0 

0 

Ref. Chord (m) 

4 

Ref. semi-span (m) 

3 

Ref. Area (m 2 ) 

7.5 

X is the quarter-chord sweep angle 


Folding Wing 


Initial 

Final [ 

Panel # 

1 

2 

1 

2 

Cr(meter) 

4 

1 

4 

1 

Ct(meter) 

1 

1 

1 

1 

AR 

1.2 

3 

1.2 

3 

X (deg) 

37 

0 

37 

0 

Fold Angle (deg) 

0 

0 

0 

180 

Ref. Chord (m) 

4 

Ref. semi-span (m) 

6 

Ref. Area (m 2 ) 

10.5 


Table 1. Telescoping wing parameters. Table 2. Folding wing parameters. 


Bat Wing 


Initial 

Final ! 

Panel # 

1 

2 

3 

4 

1 

2 

3 

4 

Cr(meter) 

0.66 

0.4 

0.4 

0.13 

2.77 

1.45 

1.45 

0.13 

Ct(meter) 

0.4 

0.4 

0.13 

0.13 

1.45 

1.45 

0.13 

0.13 

AR 

1.42 

1.89 

2.84 

5.67 

0.18 

0.26 

0.48 

2.89 

X (deg) 

5.04 

10 

5.04 

10 

40.9 

60 

40.9 

60 

Fold Angle (deg) 

0 

0 

0 

0 

0 

0 

0 

0 

Ref. Chord (m) 

0.66 

Ref. semi-span (m) 

3 

Ref. Area (m 2 ) 

1.986 


Table 3. Bat-like wing parameters. 


Mass, kg Ixx, kg-m 2 Iyy, kg-m 2 Izz, kg-m 2 


Telescoping and Folding 
Wings 

Center Body 

985.13 

1191.2 

3176.1 

2119.5 

Outer Wing 

58.7 

44.5 

0.75 

44.5 

Bat-like Wing (10° deg Configuration) 

985.13 

1191.2 

3176.1 

2119.5 


Table 4. Mass and inertia properties of telescoping, folding, and bat-like wings. 
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IV. Data Transfer among Disciplines 

Accurate and automated data transfer among various disciplines is one of the key ingredients for the loosely 
coupled analysis and design process. Often discipline models share the same geometry, but the models (e.g., 
aerodynamics and structures) have dissimilar meshes. In addition, the data transfer process between disciplines may 
be subjected to additional constraints, such as conservation of forces, moments, and energy. One source of error 
often occurs when models have dissimilar levels of geometry details. For example, an aerodynamic mesh generally 
resembles the true geometry of the aircraft; the mesh includes details such as pylons, nacelles, flaps, and slats. 
However, a structural mesh generally represents only the major structural modules, such as a wing box, but flaps and 
slats either are represented as simple beam elements or are completely excluded. The work by Samareh and Bhatia 10 
provides a general algorithm for data transfer that is accurate and satisfies the required constraints. Figure 6 shows 
the results of load transfer for a morphing configuration: (a) aerodynamic mesh, (b) coefficients of pressure on the 
aerodynamic mesh, (c) load vectors on the aerodynamic mesh, (d) a simple mesh for multibody analysis model 
(colors corresponds to different wing panels), and (e) transferred load vectors on the multibody mesh. The load 
transfer process used in this study is accurate and guarantees conservation of forces and moments. 



Fig. 6 Load transfer for a morphing vehicle. 


V. Aerodynamic Database Generation 

One key assumption in creating the aerodynamic database is that the shape changes occur at much slower rate 
than the fluid particles traveling across the wing chord. This assumption allows one to use a series of steady 
aerodynamic solutions to model the flow field. For this, a MATLAB based software coupled with the potential flow 
aerodynamic code CMARC is used. MATLAB is a high performance software application for engineering and 
scientific computation and visualization, and CMARC code is a commercial version of the potential-flow panel 
method code, PMARC 12 , developed at NASA Ames. The CMARC code is capable of calculating the surface 
pressure distribution around a 3 -dimensional body, assuming that the flow field is incompressible, inviscid, and 
irrotational. It also has the capability to calculate skin friction over the wing; however, this capability was not used 
for this study. This section presents the process of aerodynamic mesh generation, CMARC execution, and 
aerodynamic database generation. The process can be easily extended for aerodynamic analysis based on 
computational fluid dynamics (CFD). 
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The aerodynamic database used in this study has as input variables the angle-of-attack and geometric parameters 
describing the morphing shape changes. These geometric parameters controlling the shape changes are the retraction 
distance for the telescoping wing, a folding angle for the folding wing, or a leading-edge sweep angle for the bat-like 
wing. Outputs from the aerodynamic database consist of forces and moments acting at discrete points on the wing. It 
is at these discrete points that the multibody dynamic analysis model (mesh) accepts external loads. For example, for 
DYMORE models, the discrete points in the database are located along the wing quarter-chord line. For LMS 
Virtual Lab (VL), the discrete points are located on the leading and the trailing edges of the wing. 

Figure 7 shows a planform view of a bat-like in-plane morphing wing similar to concepts developed by 
NextGen Aeronautics and a sample input script required by MATLAB. This input script describes a four-section 
bat-like wing geometry with corresponding root chord, tip chord, aspect ratios, quarter sweep angles, and fold angles 
between each section. For the example shown in Fig. 7, each wing section has 10 panels in the spanwise direction. 
There are 50 points (25 points on the top and 25 on the bottom) per airfoil section and the wing configuration option 
allows for a full- wing geometry or a half-wing geometry to be executed in CM ARC. Based on this input, the 
MATLAB script generates a surface mesh in the format required by CMARC. This MATLAB mesh generation 
script offers great flexibility to change any or all the geometric parameters and to generate CMARC solutions 
quickly, without user interaction, for a vehicle with large shape changes. 

#WING CREATION INPUT DECK 
^number of sections 
4 

#root chord lengths 
0.66 0.4 0.4 0.13 
#tip chord lengths 
0.4 0.4 0.13 0.13 
#Aspect Ratios 
1.42 1.89 2.84 5.67 
^Quarter Sweep Angles (degrees) 

5.04 10.0 5.04 10.0 

^Number of panels comprising each wing 
10 10 10 10 

#Fold angles (degrees) 

0.0 0.0 0.0 0.0 
#Spline points 
25 

#Wing configuration option (1 full, 2 half) 

1 

#Solver type option (1 CMARC, 2 Vorview) 

1 

#Units option (2 for English (feet), 1 for Metric) 

1 

#Mach number (For flow solver) 

0.2 

#Alpha schedule vector (For flow solver) 

-10 -5 0 5 10 

^Reynolds number based on chord and dimensional kinematic viscosity 
480000 0.0001 

#File consisting of airfoil point definition 
naea001,2 . dat 









dr 

T "1 

r r 

dr 
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3 
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1 2 
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Fig. 7 Bat-like wing planform and sample input script for 
MATLAB script to generate bat wing geometry. 


http://www.nextgenaero.com 
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In addition to CM ARC, the MATLAB script has a link to a vortex lattice code called Vorview 13 (an option 
available in the input script under solver type option). The MATLAB script reads the variables describing the initial 
and the final morphing shapes from the geometry input deck and linearly interpolates for the intermediate shapes 
from the initial and the final values of leading edge sweep angles. Then, MATLAB creates a CM ARC mesh and 
invokes CMARC to obtain the aerodynamic solution. Figure 8 shows sample mesh and surface pressure solution 
results for the bat-like wing shown in Fig. 7. It took one second to generate each mesh and 1 1 seconds to generate 
each CMARC solution on a 3.2 GHz PC running Linux OS. 



Fig. 8 Meshes and coefficients of pressure for a bat-like wing, 
(left: un-morphed, middle: 50% morphed, right: fully-morphed). 


VI. Multibody Dynamics Analysis 

Classical finite element-based structural analysis methods have proven to be robust design tools for the 
determination of static, dynamic, aeroelastic, and aeroservoelastic response of conventional aircraft configurations. 
However, the structural and aeroelastic analysis of morphing vehicles capable of radical shape changes require tools 
that enable the modeling of integrated mechanisms, structures, controls, and aerodynamic forces. Also important is 
the ability to analyze a complete virtual aero-structural morphing system across its full dynamic range of motion 
with large geometric and structural deformations without having to develop time consuming finite element models 
of a succession of intermediate morphed configurations. Existing multibody dynamic analyses offer the ability to 
model integrated aero-structural-control characteristics of morphing vehicles in the time domain when coupled with 
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an aerodynamic model to aid in rapid prototyping and design optimization of various morphing concepts. This paper 
focuses on the use of two multibody dynamic analysis codes known as LMS Virtual Lab (VL) and DYMORE. 

A. LMS Virtual Lab and MATLAB Simulink Implementation 

To create the VL model for both folding and telescoping configurations, the outer mold line is read in from the 
computer-aided design (CAD) geometry, which was generated by the geometry process. Each moving surface, like 
the right and left wing, is modeled as an independent body connected using joints. For a telescoping wing concept, a 
translational joint is used to connect each wing-section to the center body. For the case where the wings fold inward, 
a revolute joint is used instead. To incorporate external aerodynamic loads into the model, eighteen points are 
selected on the vehicle surface where the loads are applied, 6 on the center body and 6 on each wing panel. These 
aerodynamic loads act in all three directions at each point. Coupling of the multibody code with the aerodynamic 
loads is accomplished within the MATLAB/Simulink software environment and the co-simulation capabilities of 
VL. The block diagram, shown in Figure 9, shows the model elements interconnections. To simulate the vehicle VL 
creates the multibody equations of motion with inputs in terms of aerodynamic forces, control surface commands, 
and support constraints, while MATLAB couples and integrates the equations in time to compute vehicle response. 
In the configuration studied, the vehicle is mounted on a pedestal free to rotate. During simulations, the desired 
angle of attack is reached by controlling reaction torques from the vehicle onto the mounting pedestal. Although this 
condition is not consistent with free flight, it is consistent with many wind tunnel testing. The aerodynamic database, 
on the other hand, accepts control surface commands (geometric parameters) and angle of attack values to produce 
the aerodynamic loads that are feedback to the multibody model. Although it is possible to loosely couple models 
like this at each flight condition manually, the planform changes are so radical that this is the most efficient way to 
study these types of vehicles. 




Fig. 9 Simulink/Virtual Lab block diagram of telescoping and folding wing models 
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B. DYMORE Implementation 

DYMORE 11 , developed by Professor Olivier Bauchau of Georgia Institute of Technology, is a finite element- 
based tool for the analysis of non-linear elastic multibody systems with arbitrary system topologies. The DYMORE 
element library consists of beams, shells, rigid bodies, joint constraints, contact joints, backlash/freeplay joints, and 
actuator elements. Sensors can be assigned to various elements to observe deformations, displacements, forces, and 
moments. Although not used in this study, DYMORE includes a state-space 2-D unsteady lifting line aerodynamic 
model and the ability to tightly couple the structural model to the Navier-Stokes structured grid flow solver 
OVERFLOW via structured files. Aeroservoelastic problems can also be formulated within the DYMORE 
environment through controller elements, which couple model sensors to joint/actuator motions via linear and user- 
defined control law definitions. Thus, morphing mechanisms of varying degrees of complexity can be modeled 
using DYMORE with the inclusion of elastic, aeroelastic, and aeroservoelastic affects 14 . 

For the present study, DYMORE models of generic telescoping, folding, and bat-like in-plane morphing vehicle 
concepts have been formulated and coupled to the MATLAB steady aerodynamic database described in the previous 
section. Using DYMORE’ s Fluid Structure Interface (FSI) formatted file exchange architecture, commanded 
morphing displacements, rotations, MATLAB/CMARC steady forces, and moments at the quarter chord are 
exchanged at each time step. The structure and morphing mechanisms of the three morphing concepts described in 
this paper are modeled in DYMORE as beam finite elements connected by translational prismatic joints and/or 
rotational revolute joints. The present DYMORE models are restrained in all translational degrees of freedom, pitch, 
and yaw, but are allowed to rotate in roll about a re volute joint to explore the use of asymmetric morphing actuation 
for roll control. Although structural deflections, rotations, and velocities of beam element nodes are passed to the 
MATLAB aerodynamic database at each time step, these are not utilized in the steady aerodynamic database of this 
study but could be utilized by an unsteady aerodynamic code as part of an aeroelastic or aeroservoelastic analysis. 

VII. Discussion of Simulation Results 

This section presents the results for three wings undergoing large shape changes: telescoping wing, folding 
wing, and bat-like wing as shown in Fig. 3. The aerodynamic results for these three cases are presented in the 
Appendices A-C. For each of the three cases, the entire setup and the analysis process took approximately 3-5 days. 

A. Telescoping Wing 

To compare results from VL and DYMORE when undergoing identical shape changes, a change schedule was 
developed and employed within both multibody codes to produce asymmetric planform changes with 1-cos 
telescoping schedule. Figure 10 shows a schematic with the orientation of the joint force F z , moment M x , and 
retraction 8. Results for the telescoping concept are shown in Fig. 11; the top plot is the right wing commanded 
retraction that ranged from 0 to 3 meters, solid blue is VL and dashed green is DYMORE. All results are obtained 
with an angle of attack of 2.47 degrees, dynamic pressure Q=2837 Pa, Mach=0.2, and lg conditions. In order to 
induce vehicle roll, the right and left wing retraction 
commands are slightly out of phase and the 
difference is shown as “A=5 r -5” underneath the top 
of plot. This command differential produced a 
maximum vehicle roll of 46 degrees in both 
simulations with a 10-12 % difference between 
solutions. Similarly, the predicted maximum roll 
rate is 3.6 deg/sec with a 10-14% difference 
between solutions. Another important aspect of 
aircraft undergoing large shape changes is joint 
loading. For the telescoping case, the bottom two 
plots in Fig. 11 show the normalized vertical shear 
force F z l(QA ref ) and torque M x l(QL ref A ref ) at the 

right wing joint measured in local coordinate axes. 

The reference quantities are listed in Tables 1 and 4. 

When evaluating the loads, differences in the 
aerodynamic discretization approaches between VL 
and DYMORE are more apparent. For example, VL 
uses 6 points on each outer wing to apply the 

Fig. 10 Free body diagram for VL and DYMORE. 
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aerodynamics loads versus 1 1 in DYMORE. Nonetheless, the net forces and moments about the center body are the 
same. During telescoping the force amplitudes are modulated as the exposed wetted area is retracted inside the 
center body. This process introduces step discontinuities in the force and moments that are more apparent in the VL 
because fewer points are used. Furthermore, because of the load discretization in VL wing forces at the junction 
between the right/left wing and the center body are all applied to the body side (as opposed to the wing), resulting in 
a step discontinuity at time zero when compared to DYMORE results. This discontinuity has little or no effect on 
roll and roll rate. 



Fig. 11 Comparison of telescoping wing results 
using VL and DYMORE. 
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B. Folding Wing 

Results for a folding wing concept are obtained using a shape change schedule similar to that used in the 
telescoping case 1-cos function, but now the left and right wing are commanded to rotate 90 degrees in 34 seconds. 
At the top of Fig. 12 is the commanded rotation schedule in degrees for the right wing and under it is the difference 
A between the right and left rotation commands (A now used to refer to angular differences) showing a maximum 
difference of 1.5 degrees. Note that even this small difference causes a maximum roll of 18.8 degrees roll in 
DYMORE and 17.6 degrees in VL with corresponding maximum roll rates of 10.3 deg/sec and 8.65 deg/sec, 
respectively. Although the results from both modeling approaches are close, it is apparent that aerodynamic load 
discretization is the main cause of the slightly different values. When comparing normalized shear forces and 
actuation moments at the joint, the forces are also within the aerodynamic database load discretization error between 
VL and DYMORE. In addition, the spikes in the VL load plots are due to coupling of the high-bandwidth pedestal 
support control system used to maintain the pitch and yaw attitude during roll simulation. 
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Fig. 12 Comparison of folding wing results using DYMORE and VL. 
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C. Bat-Like Wing 

A DYMORE multibody model of the bat-like wing in-plane morphing concept has been constructed and coupled 
to the bat-like wing aerodynamic database discussed in the previous section. Figure 13 shows the DYMORE 
multibody model. It also shows the connectivity of beam finite elements to revolute and prismatic sliding joints to 
allow the scissor-like in-plane morphing motion, which greatly affects both sweep and wing planform area. The 
right and left sides of the model can morph independently to allow for the investigation of asymmetric morphing 
commands for roll control. Results are presented in Fig. 14 for an angle of attack of 2.47 degrees, dynamic pressure 
Q=2837 Pa, Mach=0.2, and lg conditions. Similar to the telescoping and folding wings, an asymmetric 1-cos 
schedule of leading-edge sweep angle from 10 to 60 degrees is applied as shown in the top two plots of Fig. 14, 
which illustrate the right wing sweep and resulting “A” (right wing minus left wing sweep angles). 



Fig. 13 DYMORE multibody model of the Bat-like wing in-plane morphing concept. 


Figure 14 also illustrates the resulting roll 
motion due to asymmetric morphing sweep 
schedule. Right and left linear actuation forces 
that are required to maintain right and left wing 
morph positions while under gravitational and 
aerodynamic forces are shown at the bottom of 
Fig. 14. These actuation forces arise due to the 
component of aerodynamic lift in the plane of 
the wing at angle-of-attack as illustrated in Fig. 
15. Additional actuation forces would be 
required to overcome friction in all joints in 
vehicle. However, these frictional joint forces 
are not currently modeled in DYMORE. 




Fig. 15 In-plane and out-of- Fig. 14 DYMORE results for the Batwing concept, 

plane components of lift. 
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VIII. Concluding Remarks 

This paper introduced a simple and efficient integrated approach to conduct multidisciplinary analysis for 
aircraft with large shape changes. The integrated approach has reduced the set up time from months to days, and the 
approach can handle diverse morphing concepts. The analysis modules are loosely integrated to allow for inclusion 
of existing off-the-shelf analysis modules, and this feature was demonstrated with the use of two different multibody 
dynamics software. Results from the two multibody dynamic analyses software show good correlation between the 
vehicle roll, roll rates, and critical load time histories at the joint for telescoping and folding morphing vehicle 
concepts. Although the two multibody simulation programs were significantly different, roll and roll rates were 
predicted within 10-12%. 

As a result of this study, it is clear that careful attention must be given to the distribution of aero-structural load 
interfaces to minimize the effects of discontinuities in aerodynamic load application for some morphing 
configurations. On the other hand, gross vehicle performance predictions were in good agreement even with slightly 
different load distributions and update rates in the multibody codes. Also demonstrated was the use of parametric 
models to simplify and automate the model generation process for rapid evaluation of three diverse morphing 
concepts within the same environment. Our approach needs to be modified for cases that need to include model 
flexibility and time-accurate aerodynamics. 
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Appendix A: Aerodynamic Results for the Telescoping Wing 

Telescoping the outer wing section has two major benefits; it reduces the wing loading, and the variable aspect 
ratio improves the aerodynamic performance. By adjusting the amount of telescoping, one can control the lift and 
drag characteristics for various missions. The planform for telescoping model used here is shown in Fig. 4, and the 
corresponding parameters are listed in Tables 1 and 4. All reference quantities are based on the un-morphed shape. 
Figure A1 shows the meshes and CM ARC solutions for the telescoping wing at different stages of morphing. As 
described in section V, MATLAB is used to interpolate the aspect ratio of the outer wing section from its initial and 
final values. Then, MATLAB 
creates the aerodynamic 
database as a function of the 
angle-of-attack and the 
telescoping distance. Fifty 
five telescoping CMARC 
models (11 shape changes 
and 5 angles of attack) are 
computed and stored in the 
aerodynamic database. 

Figure A2 shows contour 
plots for the integrated 
coefficients for lift (C L ), 
pressure drag (C di ) and 
moment (C m ) as functions of 
angle of attack and 
telescoping distance. 

Figure A3 shows the 
aerodynamic forces acting on 
the initial and the final configurations of the telescoping wing as used in VL and DYMORE analyses. The red color 



Fig. A1 Meshes and coefficients of pressure for a telescoping wing, 
(left: un-morphed, middle: 50% morphed, right: fully-morphed). 


CL 


Cdi 


Cm 





Fig. A2 Lift, pressure drag, and moment coefficients for a telescoping wing. 


corresponds to forces at the center body section (stationary), while the blue color corresponds to the outer wing 
section (morphing section). At the top of the figure is the un-morphed configuration and the fully-morphed 
configuration is shown at the bottom. For the un-morphed configuration, the outer telescoping section of the wing is 
inside the center body and has no forces acting on it. For the final morphed configuration, forces act on the center 
body and the outer wing. For DYMORE, the forces are applied to the points located along the wing quarter-chord 
line. Figure A4 shows the forces used in VL at various angles of attack and amount of telescoping. The red color 
corresponds to forces at the center body section (stationary), while blue corresponds to the outer wing section 
(morphing section). Because the structural elements in the center body section are larger, the nodal points on the 
center body section have larger nodal forces. 
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a) Forces for DYMORE analysis. 



Outer wtng 


Outer wing 


b) Forces for Virtual Lab analysis. 


Fig. A3 Initial and final forces for multibody dynamic analyses of a telescoping wing. 
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Appendix B: Aerodynamic Results for the Folding Wing 

Folding the wing outer section, as shown in Fig. Bl, reduces wing area and increases the effective sweep angle. 
These changes— area and sweep angle— create a vehicle that is optimum for two distinct missions. The planform for 
the folding model used here is shown in Fig. 5, and the corresponding parameters are listed in Tables 2 and 4. As 
before, all reference quantities are based on the un-morphed shape. Figure Bl shows the meshes and CMARC 
solutions for the folding wing at different stages of morphing. Similar to the telescoping case, MATLAB 
interpolates for the fold angle from its initial and final values and creates the aerodynamic database as a function of 
the angle-of-attack and the fold angle. Figure B2 shows the integrated coefficients for lift C L , pressure drag Cdi, and 



Fig. Bl Meshes and coefficients of pressure for a fold wing, 
(left: un-morphed, middle: 50% morphed, right: fully-morphed). 



Fig. B2 Lift, pressure drag, and moment coefficients for a folding wing. 


moment C m as functions of angle of attack and fold angle. Figure B3 shows forces acting on VL model at various 
angles of attack and folding angles. 
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Fig. B3 Forces acting on folding configuration for VL. 


Appendix C: Aerodynamics Results for the Bat-Like Wing 

Bat-like wing allows area and sweep angle changes that are optimum for two distinct missions. The planform for 
bat-like model used here is 
shown in Fig. 4, and the 
parameters are listed in Tables 3 
and 4. We used the initial root 
chord, the wing area, and the 
semi span for the aerodynamic 
reference quantities. Again, all 
reference quantities are based 
on the un-morphed shape. 

Figure Cl shows the meshes 
and CMARC solutions for the 
bat-like wing at different stages 
of morphing. MATLAB 
interpolates the wing planform 
parameters from its initial 
values to its final values of the 
leading edge sweep angle and 
creates an aerodynamic 
database, which is a function of 
the angle-of-attack and the 
sweep angles. 



Fig. Cl Meshes and coefficients of pressure for a bat-like wing, 
(left: un-morphed, middle: 50% morphed, right: fully-morphed). 
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Figure C2 shows the integrated coefficients for lift C L , pressure drag C di , and moment C m as functions of angle of 
attack sweep angle. Figure C3 shows forces acting on DYMORE model at various angles of attack and bat angle. 
The figure also highlights the beam element nodes and aerodynamic airstation points where aerodynamic forces and 
moments from MATLAB database are applied. 


CL Cdi Cm 



Fig. C2 Lift, pressure drag, and moment coefficients for a bat-like wing. 
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Fig. C3 Forces acting on a bat-like wing configuration for DYMORE model. 
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